Similarity theory and calculation of turbulent fluxes at the surface for the stably stratified atmospheric boundary layers

In this paper we revise the similarity theory for the stably stratified atmospheric boundary layer (ABL), formulate analytical approximations for the wind velocity and potential temperature profiles over the entire ABL, validate them against large-eddy simulation and observational data, and develop an improved surface flux calculation technique for use in operational models.Comment: The submission to a special issue of the Boundary-Layer Meteorology devoted to the NATO advanced research workshop Atmospheric Boundary Layers: Modelling and Applications for Environmental Securit

Global Intermittency and Collapsing Turbulence in the Stratified Planetary Boundary Layer

Direct numerical simulation of the turbulent Ekman layer over a smooth wall is used to investigate bulk properties of a planetary boundary layer under stable stratification. Our simplified configuration depends on two non-dimensional parameters: a Richardson number characterizing the stratification and a Reynolds number characterizing the turbulence scale separation. This simplified configuration is sufficient to reproduce global intermittency, a turbulence collapse, and the decoupling of the surface from the outer region of the boundary layer. Global intermittency appears even in the absence of local perturbations at the surface; the only requirement is that large-scale structures several times wider than the boundary-layer height have enough space to develop. Analysis of the mean velocity, turbulence kinetic energy, and external intermittency is used to investigate the large-scale structures and corresponding differences between stably stratified Ekman flow and channel flow. Both configurations show a similar transition to the turbulence collapse, overshoot of turbulence kinetic energy, and spectral properties. Differences in the outer region resulting from the rotation of the system lead, however, to the generation of enstrophy in the non-turbulent patches of the Ekman flow. The coefficient of the stability correction function from Monin-Obukhov similarity theory is estimated as (Formula presented.) in agreement with atmospheric observations, theoretical considerations, and results from stably stratified channel flows. Our results demonstrate the applicability of this set-up to atmospheric problems despite the intermediate Reynolds number achieved in our simulations. © 2014 The Author(s)

Evaluation of urban local-scale aerodynamic parameters: implications for the vertical profile of wind speed and for source areas

Nine methods to determine local-scale aerodynamic roughness length (z0) and zero-plane displacement (zd) are compared at three sites (within 60 m of each other) in London, UK. Methods include three anemometric (single-level high frequency observations), six morphometric (surface geometry) and one reference-based approach (look-up tables). A footprint model is used with the morphometric methods in an iterative procedure. The results are insensitive to the initial zd and z0 estimates. Across the three sites, zd varies between 5 – 45 m depending upon the method used. Morphometric methods that incorporate roughness-element height variability agree better with anemometric methods, indicating zd is consistently greater than the local mean building height. Depending upon method and wind direction, z0 varies between 0.1 and 5 m with morphometric z0 consistently being 2 – 3 m larger than the anemometric z0. No morphometric method consistently resembles the anemometric methods. Wind-speed profiles observed with Doppler lidar provide additional data with which to assess the methods. Locally determined roughness parameters are used to extrapolate wind-speed profiles to a height roughly 200 m above the canopy. Wind-speed profiles extrapolated based on morphometric methods that account for roughness-element height variability are most similar to observations. The extent of the modelled source area for measurements varies by up to a factor of three, depending upon the morphometric method used to determine zd and z0

The Critical Richardson Number and Limits of Applicability of Local Similarity Theory in the Stable Boundary Layer

Measurements of atmospheric turbulence made over the Arctic pack ice during the Surface Heat Budget of the Arctic Ocean experiment (SHEBA) are used to determine the limits of applicability of Monin-Obukhov similarity theory (in the local scaling formulation) in the stable atmospheric boundary layer. Based on the spectral analysis of wind velocity and air temperature fluctuations, it is shown that, when both of the gradient Richardson number, Ri, and the flux Richardson number, Rf, exceed a 'critical value' of about 0.20 - 0.25, the inertial subrange associated with the Richardson-Kolmogorov cascade dies out and vertical turbulent fluxes become small. Some small-scale turbulence survives even in this supercritical regime, but this is non-Kolmogorov turbulence, and it decays rapidly with further increasing stability. Similarity theory is based on the turbulent fluxes in the high-frequency part of the spectra that are associated with energy-containing/flux-carrying eddies. Spectral densities in this high-frequency band diminish as the Richardson-Kolmogorov energy cascade weakens; therefore, the applicability of local Monin-Obukhov similarity theory in stable conditions is limited by the inequalities Ri < Ri_cr and Rf < Rf_cr. However, it is found that Rf_cr = 0.20 - 0.25 is a primary threshold for applicability. Applying this prerequisite shows that the data follow classical Monin-Obukhov local z-less predictions after the irrelevant cases (turbulence without the Richardson-Kolmogorov cascade) have been filtered out.Comment: Boundary-Layer Meteorology (Manuscript submitted: 16 February 2012; Accepted: 10 September 2012